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^ ■ Abstract 
00 

\ When relating the strong coupling as, measured at the scale of the Z boson 

' mass, to its numerical value at some higher energy, for example the scale of Grand 

■ Unification, it is important to include higher order corrections both in the running of 

Xf^ \ Us and the decoupling of the heavy particles. We compute the two- loop matching 

' coefficients for as within the Minimal Supersymmetric Standard Model (MSSM) 

which are necessary for a consistent three-loop evolution of the strong coupling 
Qh, constant. Different scenarios for the hierarchy of the supersymmetric scales are 

considered and the numerical effects are discussed. We find that the three-loop 
effects can be as large as and sometimes even larger than the uncertainty induced 
by the current experimental accuracy of as(M^). 
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1 Introduction 



In spite of describing all current experimental precision data extremely well, it is com- 
monly believed that the Standard Model (SM) is not the ultimate theory of particle 
physics. This is due to deficiencies like the so-called fine tuning problem which arises 
from quadratic divergences combined with the large difference between the electro-weak 
and the Planck scale. 

A popular extension of the SM that cures this problem is the so-called Minimal Su- 
persymmetric Standard Model (MSSM). In addition, one of its most attractive features 
is the possibility of unification of the strong, the weak, and the electro-magnetic coupling 
constants. While having very different numerical values at the electro-weak scale, they 
may have a common intersection point at a scale /igut ~ 10^^ GeV [1-6] which is there- 
fore considered to be the typical scale of a Grand Unified Theory (GUT). In the minimal 
SM, such an intersection can not be achieved. 

The relation of the coupling constants at different energy scales is governed by the 
renormalization group running, accompanied by the appropriate decoupling of heavy par- 
ticles, to be described in more detail in Section |21 The most recent studies of the coupling 
constants at the GUT scale as evolved from their experimental values at the weak scale 
have been performed in the context of the SUSY Parameter Analysis project (SPA) [8]. It 
is the purpose of this paper to determine the precise evolution of the strong coupling Ug- 
In particular, we evaluate the two-loop decoupling relations within the MSSM. In combi- 
nation with the known beta function, this allows us to determine ^^(/iGUT) at three- loop 
accuracy, and to estimate its theoretical uncertainty. 

The outline of this paper is as follows: in Sections El and 01 we review the theoret- 
ical framework and introduce our notation. In Section we quote the renormalization 
constants and the j3 function in the DR scheme within the MSSM. Our results for the 
decoupling relations are given in Section and the phenomenological implications are 
discussed in Section El Our findings are summarized in Section [71 

2 Decoupling of heavy particles 

In the framework of QCD and its supersymmetric extension, it is convenient to use a 
mass-independent renormalization scheme such as MS or DR. Thus, by construction, the 
beta function governing the running of Us is independent of the particle masses. It only 
depends on the particle content of the underlying theory, i.e., the number of active quarks, 
squarks and gluinos, denoted by n/, n^, and rig in what follows (a more precise definition 
of these labels is given subsequent to Eq. (HHI))- 

The price for this simplicity is that, in contrast to momentum subtraction schemes, 
the Appelquist-Carrazone decoupling theorem [9] does not hold. This is well-known in 
pure QCD where, at an energy scale /i. Green functions involving heavy quarks of mass 
mh depend logarithmically on ii/mh- In order to avoid potentially large logarithms for 
yU <^ ruh one has to decouple the heavy quark from the theory [10]. As a consequence, one 
obtains different coupling constants in the various energy regimes. They are connected 
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by the so-called decoupling or matching relations. This is usually indicated by labeling 
Us with a superscript which determines the number of active particles that contribute to 
its running. Considered as a function of /i, shows a kink at one-loop and even a jump 
at higher orders when changing from one energy regime to the other. 

The precise value where the heavy degrees of freedom are integrated out is not fixed by 
theory. Since the explicit dependence on this matching scale is logarithmic, it is natural 
to choose /i ~ TTih- The dependence of theoretical predictions on the variation of the 
matching scale around this intuitive value can be used as an estimate of the theoretical 
uncertainty. Reduction of this uncertainty can only be achieved by means of higher order 
calculations. 

Within QCD, the decoupling relations have been computed up to the three-loop or- 
der [11]. In the MSSM, only the one-loop relation for is known [12,13]. It is the purpose 
of this paper to extend this calculation to two loops, such that a consistent three-loop 
evolution of to the GUT scale can be performed. The matching coefficients that relate 
different energy regimes (with different values of n^, n/, and rig) involve the masses of 
the heavy particles. Whereas at one-loop order the dependence is only logarithmic, it in- 
volves much more complicated functions at two-loop order. Furthermore, it is important 
whether the heavy particles have approximately similar masses and can thus be integrated 
out simultaneously, or whether the mass differences are big and the decoupling has to be 
performed step-by-step. 

In this paper we will consider different scenarios concerning the hierarchy among the 
relevant particle masses. They are defined by the following conditions: 

(A) rriu, . . . , > m^-, m^, rrit > mt 

(B) niu, mi, mi > m^, nit > "^6 

(C) iriu, mi, mi, mg mt mb 

(D) mu, mi, mi, m-g, mt m^ , 

where m^, . . . , mi are the squark masses, mg is the gluino mass and mt, mb are the top 
and the bottom quark mass, respectively. 

In all cases it is understood that masses which are separated by a comma are of the same 
order of magnitude. Note that case (B) corresponds to so-called spht SUSY which has been 
receiving much attention [14, 15] recently. Furthermore, in case (C) all supersymmetric 
masses are considered as heavy with respect to the SM ones. Squark masses are always 
taken to be larger than the top quark mass. 

If all occurring mass scales are different, the resulting formulae are quite complicated 
and unhandy. For the purpose of this paper, we thus restrict ourselves to the simplified 
scenarios where either all masses are identified (if possible), or where an expansion in 
mt/M can be performed, with M G {m^, . . . , mi, mg}. Further details on our approxima- 
tions will be given in Sections El and El 
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3 Technicalities 



We follow the framework defined in Ref. [11, 16] (see also Ref. [17]) which reduces the 
calculation of the n-loop decoupling relation to n-loop vacuum diagrams. The bare de- 
coupling (or matching) coefficients Ci", (i are introduced as 



ccg SI ^ ccg ' \-^J 



where gs = y/Airag is the QCD gauge coupling, is the gluon field, c"" is the Faddeev- 
Popov ghost field, and Tccg is the ghost-gluon vertex. The prime generically denotes the 
quantities in the effective theory, where the heavy particles have been integrated out and 
the bare objects are marked by the superscript "0". According to Ref. [11], and 
are obtained from the gluon propagator, the ghost propagator, and the ghost-gluon 



vertex, all evaluated at vanishing external momenta, (g is then determined from the 
relation 



C° = (2) 

^9 AO /To ^ ' 

S3 V S3 



The renormalized decoupling coefficient is obtained with the help of the renormaliza- 
tion constants defined through 



9°s = f^^'ZgQs, 



Go, a / /7 /^a 

^ — V^3^^J.^ 

^0,a _ . / 7 „a 
C ~ V ' 

^ccg = ZiVccg , (3) 

where /x is the renormalization scale and D = 4 — 2e is the number of space-time dimen- 
sions. From Eq. Q one obtains a relation that determines the renormalization constant 
of the strong coupling: 

We have explicitly computed the two- loop expression for Zg by evaluating Z3, and Zi 
to the corresponding order in the MSSM. The result will be given in Section 0] 
The finite decoupling coefficient is obtained via 

C. = §C, (5) 



Z'g 
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where Z'^ corresponds to the renormahzation constant in the effective theory. 

In Ref. [11, 17] the framework described above has been apphed to the computation 
of the three-loop decouphng coefficients in QCD. The same method can be applied to 
the MSSM with the complication that more than one mass scale appears. In particular, 
one has to deal with a massive gluino and different masses for the two squarks of the 
same flavour. Furthermore, the different hierarchy structures among these masses lead 
to different expressions for the decoupling coefficients. The starting point of all scenarios 
defined in the Introduction is the energy regime where all particles of the theory are 
active (i.e., they contribute to the /? function of the strong coupling); the corresponding 
coupling is denoted by ai^"^^''. The end point is given by the regime where only the five 
lightest quark flavours are active, with the corresponding coupling denoted by Os ■ In 
case (D) a one-step decoupling from ai^"^^"* to ctt' is performed. However, in order to have 
more flexibility, in scenarios (A), (B) and (C) a two-step decoupling is defined where the 
corresponding strong couplings and decoupling coefficients are defined through 



(A) 










(B) 










(C) 


«f = 








(D) 


«f = 






(6) 



The dependence on the renormahzation scale n and — for the decoupling coefficients — 
on the masses is suppressed. The main goal of this paper is the computation of the Cif 's- 

Due to the isospin constraint between the top and bottom squark masses (see, e.g., 
Ref. [18]), scenario (A) is inconsistent. However, from the theoretical point of view it 
constitutes an appealing mass pattern. Furthermore, combined with scenario (C), it can 
be used to construct more realistic settings as will be described briefly at the end of 
Sectional As our default scenario in the numerical discussion we adopt scenario (C). 

Dimensional regularization (DREG) is the method which is most widely used in order 
to regularize the divergences of loop integrals. However, it is well-known that DREG 
violates supersymmetry. A possible compromise is provided by dimensional reduction 
(DRED) [19] which enforces the equality of fermionic and bosonic degrees of freedom 
while preserving all the calculational advantages of DREG. All the results of this paper 
are thus evaluated in the DR scheme, i.e., we used DRED accompanied by minimal 
subtraction. 

The final results are expressed in terms of the strong coupling constant as in the DR 
scheme and the on-shell mass of quarks, squarks and gluinos. The experimental input, 
as{Mz), is usually given in the MS scheme, however, such that we need to convert it 
to the DR scheme. We evaluated the corresponding relation through two loops in the 
MSSM, finding 




Figure 1: Sample diagrams contributing to Z^, and Zi with gluons [g), 
ghosts (c), top quarks (t), top squarks (t) and gluinos (g). The diagrams 
contributing to the decouphng coefficients must contain at least one heavy 
particle. 



where rif is the number of quark flavours and = 1 or distinguishes whether the 
gluino is active or not. Only the QCD-part of this equation will be used in what follows. 
It is obtained by setting rig = and agrees with Ref. [20]. With n/ = 5, the input 
af^{Mz) = 0.1 187 ±0.002 [21] g_ives a^{Mz) = 0.1198 ±0.002 when using the one-loop 
approximation of Eq. ((7j), and a^^{Mz) = 0.1200 ± 0.002 with two-loop accuracy. 
Note that, in what follows, we will often use the short hand notation 

as ^ af^ifi) . (8) 



4 Renormalization constants and function 

In this section we present the results for the DR renormalization constants. Since the 
DR scheme is mass- independent, the calculation can be reduced to massless integrals 
depending on only one external momentum. More precisely, for Z^, Z^ and Zi one has 
to consider the gluon propagator, the ghost propagator, and the gluon-ghost vertex; in 
the latter case, the external momentum is chosen to flow through the ghost line. Sample 
diagrams are shown in Fig. ^ An important check is provided by the QCD result which 
we re-derive for general gauge parameter ^, defined through the gluon propagator 



± ie 

Up to two-loop order our results read [22] 



^r(^) = , • (9) 
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' 12 ~ y + 8 



192 96^ ^^ 384 64/ 96 ^ ^ ^' 
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128 64 



256 128 



{uf + 2ns - nsUg] 



CaT 
64 



(lOrij — — SnsTig) 



1 + 



+ 



+- 



1 + 



+ 



TT 

TT 
1 

e 

vr 

vr 



( + 



16 64 ^ 



^ + + CAT(2nf + ris) — 

256 512/ ^ ^ ^ '^64 



43 



^ ' 768 384"^ 512 



5 e 



C^T(10n/ + lln 



384 



7e , e 



e2'"^ll28 128 64 



3 7^ C 

+ 



128 256 256 



(10) 



where nj and are the numbers of active quark and squark flavours, respectively, and 
rZg = or 1 distinguishes whether the gluino is active or not. For example, in the full 
f = Hg = 6 and Ug = 1. The terms involving ngUg stem from diagrams 



MSSM it is n 

containing a quark-squark-gluino vertex.^ The precise meaning of ag depends on the 



actual values of these constants. E.g., for nj = 5, Ug = rig - 
group SU(A^c), the colour factors in Eq. (fTUj) are given by 



it is a, 



a 



(5) 



For a gauge 



2Nr 



T = 1/2 



The combination of Z^, and Zi according to Eq. Q leads to the renormalization 



^Note that in this paper ql, qr are either both active or both integrated out. Also, we always assume 
that, if q is active, so is q {q € {u, d, c, s, t, b}). 
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constant for the strong coupling 



1 + 



TT e 



^ , 11 1 ^ T.^ 



+- 



^2 



17 



121 
144 



-n,J+-(2n,+ 



;i2) 



which is needed for the renormahzation of the decouphng coefficient (cf. Eq. ©). 
From Zg one obtains the /3 function defined through 



/i 



d as 

d/i^ TT 



-) Eftf- 

71 J ^-^ V TT 



j>0 



(13) 



with the result 



1 
4 
1 

16 



Ca 



11 2 \ 2 , 

y ~ 3^V " 3 ' ^ 



34 16 



—n-r. 



(14) 



-CAT{5nf + Ug — SusTig) — ACpTijif + 2^^ — n^ra^) 



Note that due to the terms oc UsTig this result can also be used for scenarios where the 
gluino is decoupled but (some) squarks are not. However, for the three-loop coefficient 
(32 (see below) such terms are not separately available. Thus in the following we identify 

The results for /3o and jSi are renormahzation scheme independent. The pure QCD 
part agrees with the well-known MS result which can be found in Ref. [23], for example. 
The contribution from the SUSY particles agrees with Ref. [24]. 

The two-loop decoupling relations presented in the next section have to be used in 
conjunction with the three-loop coefficient of the /? function. Subtracting the pure QCD 
result in the DR scheme [20] from the fully supersymmetric expression [24], we can distin- 
guish the contributions from squarks and gluinos and multiply them by the appropriate 
labels Hg and n^: 



P2 



1 

64 



—CaCfT 



3115 
209 



899 
27 



1439 



nf\+CAT' 



50 
27 



158 
27 ^ 



2 I 2 

+ -r^n 



9 



-Us + 



259 



nf]+ ClT {lAns + 2nf) + CpT^ 



148 



6 



-n] 
9 ^ 



(15) 
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Note that since Ref. [24] assumed Uf = Us, it is not possible to distinguish the term 
proportional to ngUf from the term. For this reason, Eq. p5|) can only be used for 
Uf = Hs or Hs = 0] unfortunately, this prevents us from discussing scenario (A) in the 
numerical applications. As expected, the pure QCD result in DR, obtained by setting 
Us = rig = in Eq. is different from the expression in the MS scheme; only the 
coefficients of the colour structures C^iT^n^ and CpTuf are the same in both schemes. 

Apart from the DR renormalization constants in Eq. (llUj). we also need the one- 
loop renormalization constants for the top quark and squark mass. The explicit on-shell 
expressions are given in Eqs. (B.5) and (B.6) of Ref. [13]. The on-shell counterterm for 
the gluino mass reads (see, e.g., Ref. [25]) 



^ = l-^\-{3C^-2Tn,)+'-cJ'- + h 
m-g TT I 4e^ ^ 2 V2 2 



^9 



q i=l 

1 m. 



2 2 2 2 2 

m„ 771% — 7711 — fnz, „ 

^ (1 + L,) 1(1 + L J + 1 '-Bl-{77il m„ m,-J 

TUg Trig 771^ " 



T^-^sin2^, [5f(m|,m„m,-J -Sfi'^(m|,m„m,-J] \ , (16) 



2 THri 

q y 



where Bq'^ is defined in Eq. (B.8) of Ref. [13]. Note that the 0{e) part of the mass 
counterterms is not needed since the one-loop result depends only logarithmically on the 
masses. In Eq. (fT^ . we have used the short-hand notation 

L^=ln — , X e {q, t, g, Qi, . . .} . (17) 
7n^ 

In addition to that, the abbreviations 



^M = ln^, L^ = ln-^, Ltm = ln^, (18) 



will be used in what follows 



5 Decoupling coefficients 

As already pointed out in Section 01 the calculation of the decoupling coefficients can be 
reduced to vacuum integrals [11,16,17]. At two-loop level, such integrals can be evaluated 
fully analytically [26] for arbitrary masses. For the purpose of this paper, however, it is 
sufficient to solve them for either a single or for two largely separated mass scales. In the 
first case, the analytical result is rather short in general, while in the second case we can 
perform an asymptotic expansion in the ratio of the two mass scales and truncate the series 
after the first few terms which leads again to a fairly compact expression. Furthermore, 
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in order to shorten the expressions even more, we present the results for the SU(3) gauge 
group. 

The Feynman diagrams are generated using QGRAF [27], and subsequently translated 
by Q2E/EXP [28,29] to the notation of MATAD [30], which evaluates the vacuum integrals. 
If the integrals contain more than one mass scale, EXP performs an asymptotic expansion 
on the diagrams before passing them on to MATAD. 

Before discussing the different scenarios, we present the pure QCD result for (g in the 
DR scheme, where the top quark is considered as heavy and the remaining (rij — 1) = 5 

quarks are massless. Up to two-loop order, the decoupling coefficient relating ai"^'' and 

(nf-l) . 

as via 



reads 



with Lt = ln(/i^/m^). Note that in the MS scheme the constant term of order is —7/24 
instead of —5/24 as in Eq. (|20|) [10,11,16,31]. The transformation of Eq. (jl9|) from the 
DR to the MS scheme confirms the fermionic two-loop coefficient of Eq. (jZj) . 

Within the MSSM, the one-loop result for the decoupling coefficient relating al to 
the strong coupling where also squarks and gluinos are active can be written in the generic 
form (see, e.g., Ref. [13]) 



(cr^)^ = 1+ 



^(SUSY) 



6 12 ^ ^'^ 2 ^ I 12 ^ *^ 

s=l \ 



S = l / 



+ O (a^) , (21) 



with the notation as defined in Eqs. ^ and (HH)), and C(2) = vrVe = 1.64493. . .. For 
the sake of consistency, also the term of order e is given; it is needed in intermediate steps 
of the calculation. The sum runs over all active squark flavours. We introduced a generic 
label "SUSY" which has to be adapted to the scenarios (A)-(D) discussed above. Note 
that due to the different mass scales involved in Eq. ()21|) . exhibits a jump even at 
one-loop order if all heavy particles are integrated out at the same scale. 

Let us now discuss the various cases deflned in the Introduction and in Eq. (jU)). The 
common large mass will be denoted by M and the common smaller mass by fh. In some 
cases, the SUSY masses will be taken to be (much) larger than the top quark mass. 
However, it turns out that the expansion in the ratio converges very fast — even for equal 
masses. Since in Section IHl the inverse decoupling coefficients are needed, we will give 
explicit results for them. The analytic expressions for C,g itself can be obtained easily by 
recursive substitution of the one-loop relations. 



10 



5.1 Scenario (A) 



mi 



M ^ 171 = mi = mg , mt 



a 



(C 



Al^|2 ^(fuU) 



a 



(5) 



(C 



is a function of M only. In fact, following the method of Ref. [11, 16, 17], the masses 
of the gluino, the top squarks and all quarks have to be set to zero in the corresponding 
Feynman diagrams,^ cf. Fig.^ The result is 



1+ 



a 



7T 12 



a 



TT 



13 r , f^s J 



48 



12 



12' 



[22] 



with Us = 5. 

The Feynman diagrams leading to do not contain the five squark flavours already 
integrated out in step Al. However, one has the top squarks, gluino and top quark as 
massive. We consider two limits: 

(A2a) m = mf. 



(.A2a)2 



(5) Q 

1 + '^h. 



+ 



IT 4 



a 



(5)' 



TT 



/ 1909 , 11^3 
I "864 



+ 



108 



17^ 89^ , 

-TT H Do H + —Lr 

12 24 V4 



(23) 



(A2b) m > mt: 



1 + 



a. 



(5) 



TT 



ht + ^L^ 



265 19 35 [1 7 
^ + 24^* + T2^'^+U * + T2^'^ 



Ttt /mtY 
36 V m / 

288 V m y 




(24) 



where 5*2 = ^Cl2(|) — 0.260434. It is interesting to note that 1/(C^^'^)^ as given in 
Eq. ()24|1 evaluated for m = mt approximates 1/(C^^^)^ to an accuracy better than 1%. 

^Note that this is no restriction on the resulting effective theory: once the heavy particles are decou- 
pled, the remaining fields can have arbitrary masses <^ M . 
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5.2 Scenario (B) 



m,-,. 



mi 



M > 



m 



a ' 



ITLt 



The analytical expression for C,^^ is identical to the one of C^^, only now we have Ug = 6. 
In the case of (^"^ there is no gluino-top-stop interaction since at this stage all squarks 
are already integrated out. As a consequence the integrations factorize into parts where 
either the gluino or the top quark is the only massive particle. One obtains 



1 + 



(5) 



TT 



-Lt + -Lg 

6 2 ^ 



+ 



a 



(5)' 



2 r 



TT 



275 19 , 

\ Lt + 3Lg 

96 24 ^ 



6 



2 ^ 



(25) 



5.3 Scenario (C) 



. . . = mi = mg = M ^ mt 



The result for is given by 



1 + 



(6) 



+ 



(6) 

as 



TT 



1 ns 
- + — 
71 \2 12 

2 



L 



M 



32 4J 



12 



2 12 



with ris = 6. C,^"^ is identical to the pure QCD result of Eq. ipU]) . 



aC2 ^ aQCD 
^9 ^9 



M 



(26) 



(27) 



5.4 Scenario (D) 



M = m^i 



. . = mi = mg , mt 

{5) = (gD)2^(fuU) 



As in the case (A2) we consider two limits: 
(Da) M = mt: 

.,(5) 



1 



1 + 



as 



TT V3 12 
91 n 



:t + t4 1 + 



(5) 



TX 



+ ' 24~12j^^^+i^ + li) 



1999 11^3 17^ 5 

1 vr H 62 Us 

864 108 12 48 



(2^ 



12 



(Db) M > mt: 




+ 



+ 



M V 48 8'^*'^y' + 36 Uf 



M 



881 ^ I \ ^'^ f^t 
~7200 ^80 / ^ 288 1 M 



5 19, 

— Us H Lt 

48 24 



+ . . . 



(29) 



where we have = 6 in both cases. We again observe very good convergence of 1/ (C^*^) 
even for M = mt] thus, for our numerical analysis in Section IHl we only take the result 
from Eq. 

At this point two remarks are in order. First, it is interesting to note that the following 
relation holds: 



C1aC2 



M2 



(30) 



which constitutes an important cross check of our calculation. Second, let us remark that 
it is possible to add the terms of involving Us to C^^'^ or C^^*^ thus take into 
account additional squarks as being integrated out together with the top squark, gluino 
and top quark. Of course, the parameter Us of has to be adjusted accordingly. 



6 Numerics 

In this section we study the phenomenological implications of the two-loop decoupling 
coefficients. In particular, we compute the value of the strong coupling Q;i^"^^''(MsusY) at 
a high scale Msusy, as evolved from the experimental input value af'\Mz). As described 
in Section 121 at higher orders in perturbation theory heavy particles have to be decoupled 
properly when crossing the corresponding particle threshold. However, the precise value 
of the scale where this is done, Hth, is not fixed by theory and can be used as a means 
to estimate the theoretical uncertainty. On general grounds one expects that, at fixed 
order perturbation theory, the relation between af\Mz) and af^"'* (Msusy) is not very 
sensitive to the choice of the matching scale as long as /xth is of the order of the heavy 
mass scale. In the following we study ai^""'* (Msusy) as a function of fith for the first three 
orders in perturbation theory. 

Let us describe the procedure for computing ai^"^^"* (Msusy) from the knowledge of 
af'\Mz) in detail for scenario (C); the other cases can be treated in complete analogy. 
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800 GeV 


0.0916 ±0.0012 


A'GUT 


0.0398 ±0.00023 



Table 1: as{fi) in the DR scheme for /i = 800 GeV and /i = /icuT assuming 
scenario (C). The starting point for the prescription described above Eq. 
is defined in Eq. (|32p . The error is due to the experimental uncertainty of 
asiMz). 



First, we compute ai (/^Sf ) from ai (Mz) by using the L-loop /? function with n/ = 5 and 
= fig = 0, and numerically solving the corresponding renormalization group equation. 

In the next step, the (L — l)-loop expression of (^'^ is used to arrive at af\iJ,f^). With 
the help of the L-loop (3 function (with Uf = 6, rig = and rig = 0) we then evaluate 
ai^\n^^) and perform the matching to as^^^\fi^^) using (^^ at (L — l)-loop order. Finally, 

ai^'^"'*(MgusY) is obtained again with the help of the L-loop /3 function after setting riy = 6, 
Ug = Q and rig = 1. We apply this procedure for L = 1, 2 and 3 and study the dependence 
on the two matching scales /i^f^^ and /i^^. Of course, in scenario (D) there is only one 
matching scale and the procedure simplifies accordingly. 

Before discussing the numerical effects, let us specify the input values used in our 
analysis. In the four scenarios defined in Section |21 we need to define a large and a small 
mass scale, M and m, and the top quark mass. We choose 

TJit = 174 GeV, 
M = 1000 GeV , 

m = 400 GeV. (31) 

For the default values of the matching scales we assume /xth,i = M and /ith,2 = 'fh. 
Furthermore, according to the discussion below Eq. ((Zj), we use 

af\Mz) = 0.120 ± 0.002 , (32) 

where the label DR has again been omitted. The GUT scale will be taken as /igut = 
lO^'' GeV. 

In order to get an impression of the numerical value for ag at various scales and the 
uncertainties induced by the input as{Mz), Tab. [T] gives ««(/!) for /i = 800 GeV and 
fi = /iGUT- These numbers are obtained from Eq. ()32|) for scenario (C) with three-loop 
running and two-loop matching. Further on, the dependence on the matching scales and 
the loop order is exemplified in Fig.|21which shows as(/iGUT) for scenario (C) as a function 
of /ith,i and /ith,25 each of them varying by a factor of ten around its default value. One can 
clearly see how the result stabilizes as soon as two-loop running and one-loop matching 
(middle plane) is included. The three-loop corrections stabilize the dependence of as on 
fith,i and /ith,2 even more (upper plane), such that the result is practically independent of 
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Q^s(/^GUt) 




40000 



Figure 2: asl/iGUx) for /igut = 10^^ GeV as a function of the two match- 
ing scales /ith,i and /ith,2- The fower (green), middle (blue) and upper (red) 
planes correspond to one-, two- and three- loop running accompanied with the 
corresponding order in the decoupling relation for the scenario (C). 



the matching scales within a wide range around their default values. Within the considered 
range of /ith,i and /ith,2, the variation of as(/iGUT) is 13.8/1.6/0.3% at l-/2-/3-loop order. 

Figs. 01 and 0] show the dependence on one of the scales when the other is kept at its 
default value. In scenario (D), where only one matching scale is present, we set all SUSY 
masses to M when the variation of /ith,i is considered (Fig. E)), and to m when /ith,2 is 
varied (Fig. EJ. A drastic improvement can be observed after including two- loop running 
and one-loop matching. When the three-loop running is included, the scale variation 
becomes almost completely flat. 

In order to study the effect of different mass configurations. Figs. El^a) and (b) compare 
the three-loop curves of Figs. El and |3 among scenarios (B)-(D) (long to short dashes). 
One observes differences of the order 0.005 which is about a factor 10 larger than the 
uncertainty originating from the current experimental uncertainty in as{Mz)- Estimating 
the theoretical uncertainty by the difference between the two- and three-loop result, we 
obtain 

as(/iGUT) = 0.0398 ±0.00023 ±0.0025 ±0.00007 . (33) 



Sas{Mz) 



masses 



th 



Another practical application of our results is shown in Fig. jUl where the strong cou- 
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Figure 3: Q;s(/iGUT) as a function of /ith,i for scenario (B)-(D). The dotted, 
dashed and full lines correspond to one-, two- and three-loop running. In 
scenario (D) we set all SUSY masses to M. 



pling at a proposed center-of-mass energy ^/s = 800 GeV of a future International Linear 
Collider (ILC) is computed. While for /ith,2 ~ 200 GeV the two- and three-loop result 
practically give the same value for a^, for /ith,2 = 2000 GeV the shift induced by the 
new three-loop terms is comparable with the uncertainty induced by the experimental 
value for as{Mz) (cf. Tab.P). This underlines once more the importance of the two-loop 
decoupling and three-loop running terms. 

Fig. [7|^a) shows l/a^ as a function of from the weak to the Planck scale, /xpL = 
10^^ GeV, where the dotted, dashed and solid lines correspond to one-, two- and three- 
loop running. It can be seen that there are kinks at 400 GeV and 1000 GeV where the slope 
of the curves changes. A closer look shows that the two- and three-loop curves actually 
are discontinuous as as jumps due to the one- and two-loop decoupling relations. In this 
plot the uncertainties from as{Mz) or the various mass patterns of the supersymmetric 
particles are hardly visible. For this reason, we show in Fig. [7|^b) the ratio as{^)/asfi{fi), 
where a<j,o(/^) is the leading order result for scenario (C) with asfi{Mz) = 0.120. The 
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Figure 4: ^^(/Ugut) as a function of /ith,2 for scenario (B)-(D). The dotted, 
dashed and full lines correspond to one-, two- and three-loop running. In 
scenario (D) we set all SUSY masses to m. 



scale /X is varied between Mz and /xpl- The dotted, dashed, and solid line correspond 
to the one-, two-, and three-loop result of scenario (C), respectively. Note that below 
Aith,2 = rh = 400 GeV the difference between the dashed and full curve is very small. 
The (green) band indicates the variation originating from the experimental uncertainty 
of Q^^{Mz) (cf. Eq. (jn2I))- The dash-dotted lines correspond to the three- loop results of 
scenario (D) where the masses are set to 2M (lower) and fh (upper), respectively. They 
thus reflect the dependence on the mass parameters of the MSSM. The corresponding 
curve for scenario (B) would lie between the two dash-dotted lines and are therefore not 
shown. 

In Fig. ISlwe combine our three-loop accuracy findings for the evolution of the inverse 
strong coupling at the GUT scale (yellow band) with the two-loop order results for the 
strong (a^, red band), the electro-magnetic (ai, blue band) and the weak coupling (a2, 
green band) (for a precise definition, see Ref. [7], and references therein). The wide error 
bands correspond to the present experimental accuracy for the gauge couplings and the 
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Figure 5: Comparison of asificvT) as a function of /ith,i (a) and /ith,2 (b) where 
for the cases (B), (C) and (D) (from long to short dashes) the three-loop results 
are plotted. 
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Figure 6: as(800 GeV) as a function of /ith,i for scenario (C). The dotted, 
dashed and full lines correspond to one-, two- and three- loop running. 
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Figure 7: (a) as(/i) as a function of /i for scenario (C). The dotted, dashed and 
full lines correspond to one-, two- and three-loop running, (b) Ratio of as(/^) 
and the corresponding leading order value corresponding to a^^{Mz) = 0.1187 
for scenario (C). The dotted and dashed lines correspond to one- and two-loop 
results, respectively. The (green) band corresponds to the uncertainty of a.^ as 
given in Eq. (jH^ . The dash-dotted lines corresponds to the three- loop result 
for scenario (D) where the SUSY masses are set to 400 GeV (upper line) and 
2 TeV (lower line). 



supersymmetric mass spectrum of the SPSla reference point. The inner bands account for 
the expected improvements in the absolute errors due to the future experimental analysis 
at LHC+LC. 

7 Conclusions 

In order to get insight into the mechanism behind Grand Unification, it is necessary to 
have precise relations among the various couplings evaluated at the electro-weak and the 
GUT scale within the expected theoretical framework, in our case the MSSM. In this paper 
we provided the two-loop decoupling coefficients which, in combination with the known 
three-loop renormalization group equation, establish the three-loop relation between the 
strong coupling at Mz and /xgut- Furthermore, they ensure the independence of the 
decoupling procedure from the scale at which the heavy particles are integrated out, up 
to higher order terms in a^. 

A large variation of ^^(/iGUT) results from the uncertainty in as{Mz)- However, as 
can be seen from Eq. ()33p and Fig. [7|^b), there is also a significant uncertainty from the 
values of the individual mass parameters. The effect of including the two-loop decoupling 
relation is a significant stabilisation of the dependence of Q!s(/iGUT) on the matching scales. 

The analysis of this paper is based on certain approximations concerning the mass 
pattern (cf. Eq. (jH)) ) within the MSSM spectrum. This results in compact formulae. It is 
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1015 10^6 
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Figure 8: Running of the inverse couplings l/oj, i = 1 (blue band), 2 (green 
band), and s (red band) with two-loop accuracy around the unification point 
defined as the meeting point of ai and a2- The wide error bands take into 
account present data and the supersymmetric mass spectrum from LHC mea- 
surements within mSUGRA. The narrow bands envisage the improvements 
expected from the LHC+LC coherent analysis [7]. The yellow band shows our 
results for l/a^ with three-loop accuracy. 



desireable to derive in a future calculation more general formulae and investigate further 
the dependence on the quark masses. 
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